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Abstract 



We study the classical and quantum theory of a class of nonlinear dif- 
ferential equations on chronology violating spacetime models in which space 
consists of finitely many discrete points. Classically, in the linear and weakly 
nonlinear regimes (for generic choices of parameters) we prove existence and 
uniqueness of solutions corresponding to initial data specified before the dis- 
chronal region; however, uniqueness (but not existence) fails in the strongly 
coupled regime. The evolution preserves the symplectic structure. 

The quantum theory is approached via the quantum initial value problem 
(QIVP); that is, by seeking operator- valued solutions to the equation of mo- 
tion with initial data representing the canonical (anti) commutation relations. 
Using normal operator ordering, we construct solutions to the QIVP for both 
Bose and Fermi statistics (again for generic choice of parameters) and prove 
that these solutions are unique. For models with two spatial points, the re- 
sulting evolution is unitary; however, for a more general model the evolution 
fails to preserve the (anti)commutation relations and is therefore nonunitary. 
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We show that this nonunitary evolution cannot be described using a super- 
scattering operator with the usual properties. 

We present numerical evidence to show that the bosonic quantum theory 
can pick out a unique classical limit for certain ranges of the coupling strength, 
even when there are many classical solutions. We show that the quantum 
theory depends strongly on the choice of operator ordering. 

In addition, we show that our results differ from those obtained using the 
"self-consistent path integral" . It follows that the path integral evolution does 
not correspond to a solution of the equation of motion. 

03.70.+k, 04.90.+e 
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I. INTRODUCTION 



Spacetimes containing closed timelike curves (CTC's) provide an intriguing environment 
for the formulation of both classical and quantum physics. Because the present is influenced 
by both the past and the future, neither existence nor uniqueness is guaranteed a priori for 
solutions to initial value problems for particles and fields on such spacetimes; these issues 
underlie many of the apparent paradoxes associated with time travel. In this paper, we 
attempt to gain insight into the initial value problem for a class of nonlinear differential 
equations (which may be regarded as toy field theories) on "spacetimes" of a type intro- 
duced by Politzer |[J. These spacetimes are defined by taking the Cartesian product of a 
number of discrete points (representing space) with the real line (time) and then imposing 
certain identifications to introduce CTC's. The simple nature of these models removes many 
technical problems and allows us to pursue the analysis to its end. 

Previous studies of classical initial value problems on chronology violating spacetimes 
have mostly focussed on linear fields [0-0] and billiard ball models 0,0^0- Deutsch flO 



has also studied examples of classical computational networks with chronology violating 
components. To the best of our knowledge there have been no studies of nonlinear fields, 
except insofar as the billiard ball models may be regarded as providing insight into the 
strongly nonlinear limit. For linear fields, it turns out that one can formulate a well posed 
initial value problem under certain conditions. Friedman and Morris |2||| have rigorously 
proved existence and partial uniqueness results for massless fields propagating on a class of 
smooth static wormhole spacetimes with data specified at past null infinity. In addition, 
they have conjectured that the initial value problem is well posed for asymptotically flat 
spacetimes with a compact nonchronal region whose past and future regions are globally 
hyperbolic whenever the problem is well posed in the geometric optics limit. It is much 
easier to prove existence and uniqueness for linear fields on certain non-smooth chronology 
violating spacetimes |yj. 

In the billiard ball case, Echeverria et al. and (for similar and more elaborate systems) 
Novikov || have shown that the initial value problem is often ill posed in the sense that the 
evolution is nonunique; moreover, Rama and Sen || have given similar examples in which 
there appears to be no global self-consistent solution for certain initial data. 

One would expect that nonlinear fields should interpolate between the behaviour exhib- 
ited by linear fields on one hand and billiard balls on the other. In this paper, we will show 
that this is indeed the case for our class of models: we prove that the initial value problem 
is well posed for arbitrary data specified before the dischronal region in both the linear and 
weakly nonlinear case, but uniqueness (though not existence) fails in the strongly nonlinear 
regime. In addition to these analytic results, we give an explicit example to demonstrate 
the lack of uniqueness for a particular value of "coupling strength". We also show that 
the evolution from the past of the dischronal region to its future preserves the symplectic 
structure. 

The loss of uniqueness for interacting systems on chronology violating spacetimes en- 
tails that classical physics loses its predictive power. Various authors have expressed the 
hope that quantum dynamics on such spacetimes might be better behaved than its clas- 
sical counterpart, with attention focussing on spacetimes possessing both initial and final 
chronal regions. Friedman et al. [fLl] considered linear quantum fields and showed that, 
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provided the classical initial value problem is well posed, the quantum evolution between 
spacelike surfaces in the initial and final chronal regions is unitary; a conclusion borne out 
by Boulware [|12|] in a Gott space example (see also |1],|5||J for related results). However, the 



situation is very different for interacting fields. Both Boulware [12] and Friedman et al. [13 



found that the S-matrix between the initial and final chronal regions fails to be unitary 
in perturbative A0 4 theory. Politzer also obtained similar perturbative results in quantum 
mechanics || and also some nonperturbative results in exactly soluble models |1J] in which 
nonunitarity also arises. It is also worth pointing out that some interacting systems do have 
unitary quantum theories 

The breakdown of unitarity raises many problems for the probability interpretation of 
quantum theory; in particular, ambiguities arise in assigning probabilities to the outcomes of 
measurements conducted before [13] or spacelike separated from [TJ|| the nonchronal region. 
There have been various reactions to these problems. Firstly, Hartle [15[] has discussed how 
nonunitary evolutions can be accommodated within the framework of generalised quantum 
mechanics, and a similar proposal has also been advanced by Friedman et al. [[13] . A second 



approach has been to "repair" the theory by modifying the evolution to yield a unitary 
theory |1|[17]]. Thirdly, Hawking [|18|] has argued that one should expect loss of quantum 
coherence in the presence of CTC's and that the evolution should be specified by means 
of a superscattering operator (i.e., a linear mapping from initial to final density matrices) 
which moreover would not factorise into a unitary S'-matrix and its adjoint. On this view, 
the quantity computed using the usual rules for the S'-matrix is not the physically relevant 
quantity and its nonunitarity is irrelevant. Rather, one should compute the matrix elements 
of the superscattering operator. Deutsch |10| has also advocated a density matrix formalism 
in the context of quantum computational networks (see also |l]]). However, this prescription 
turns out to be nonlinear in the initial density matrix JT9| ]. 

For the most part, the quantisation method employed in discussions of chronology vio- 
lation has been based on path integrals in which one sums over all consistent trajectories or 
field configurations. In this paper, we follow an operator approach based on the quantum 
initial value problem (QIVP). Namely, we seek operator valued solutions to the equation of 
motion and any consistency conditions arising from the CTC's, with initial data specified be- 
fore the dischronal region and forming a representation of the canonical (anti) commutation 
relations. If there exists a unique solution with this initial data, and the evolved data to the 
future of the dischronal region also represents the commutation relations, then we say that 
the QIVP is well posed. 

We will show that the QIVP is well posed for all linear models in our class of interest 
with both Bose and Fermi statistics. The corresponding quantum theory is unitary and 
agrees with that derived by path integral methods. In the interacting case, we prove the 
remarkable fact that (with normal ordering) the QIVP always has a unique solution and 
describe how this solution may be constructed. To obtain more specific results, we consider 
the cases in which "space" consists of either 2 or 3 points. In the 2-point model, we find that 
the unique solution to the QIVP satisfies the CCR/CAR's to the future of the nonchronal 
region (so the QIVP is well posed) and that in consequence the resulting quantum theory 
is unitary. This contrasts strongly with the corresponding path integral result (generalising 
that of Politzer |]]) in which the evolution is found to be nonunitary. In consequence, and 
because the path integral also employs normal ordering, we conclude that the self-consistent 
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path integral evolution does not generally correspond to a solution of the equation of motion. 
Given the different starting points of the two approaches this is not entirely surprising. 

In the 3-point model, we show that the QIVP is ill posed for both Bose and Fermi 
statistics because the evolved data does not satisfy the CCR/CAR's to the future of the 
dischronal region. The corresponding quantum theory is therefore not unitary. We then 
discuss the nature of this evolution in the fermionic case in order to determine whether or 
not it can be described by means of a superscattering operator. To do this it is necessary to 
translate our results from the Heisenberg picture to the Schrodinger picture. Although there 
is no unique translation prescription (as a consequence of the violation of the CAR's), we are 
nonetheless able to show that no Schrodinger picture evolution consistent with the QIVP 
solution can factorise into the product of an operator and its adjoint, lending support to 
one element of Hawking's position fI8 |. However, it also transpires that no such Schrodinger 



picture evolution can be described by a superscattering matrix as it must either increase the 
trace of density matrices or map them to non-positive operators. In this sense, the loss of 
unitarity in our model is much more radical than envisaged by Hawking. 

We also study the classical limit of our quantum theory. When this limit exists uniformly 
one recovers exactly one solution to the classical equation of motion. One might imagine 
that this limit would fail when the classical theory is nonunique; however, this is not the 
case. It appears that there are bands of "coupling strength" for which the limit does exist 
even when there are many classical solutions. We present numerical evidence to exemplify 
this behaviour including an example in which the quantum theory selects precisely one 
out of 25 classical solutions. We believe the convergence bands continue to appear as the 
coupling strength is increased indefinitely. Within the convergence bands, our quantum 
theory resolves the classical nonuniqueness; for other values of the coupling strength, it is 
arguable that no classical solution is physically relevant. 

In addition, we consider the effect of altering the operator ordering used and find that the 
solutions to the QIVP can become nonunique for large quantum numbers for non-normal 
operator ordering. We study a 1-parameter family of operator orderings for the 3-point 
model and show that the resulting quantum theories are all nonunitary. 

The paper is structured as follows. We first describe our class of chronology violating 
models in Section |lj and then study the classical initial value problem for both linear and 
nonlinear fields in Sect. [Til]. Next, in Sect. |V| we discuss the quantum initial value problem 
for our models in the absence of CTC's and demonstrate its equivalence with canonical quan- 
tisation. This serves to fix our notation and definitions for Section [V] in which we uniquely 
solve the QIVP with CTC's present, and discuss the 2- and 3-point models, showing that 
the CCR/CAR's are violated in the 3-point case. This nonunitary evolution is investigated 
in Sect. [VI] and is shown not to be described by a superscattering operator. Section |VII 
treats the classical limit, whilst Sect. |VIII| contains a brief discussion of the effect of operator 
ordering on our results. In Sect. [TX], we review the self-consistent path integral formalism, 



extending and in one instance correcting the treatment given by Politzer [l]. We use this 
formalism to compute the general (unitary) evolution for the free models, obtaining agree- 
ment with the QIVP. For the 2- and 3-point interacting models we show that the QIVP 
and path integral differ. We comment on this and other issues in the Conclusion (Sect. [X|). 
There are three Appendices. Appendix [A| reproduces our treatment of the free classical 
evolution using the methods of Goldwirth et al. B, whilst Appendix O gives a derivation 
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of the quantum evolution of the free models using the formalism of Politzer rather than 
the more direct method employed in the text. In Appendix |y, we present the details of a 
calculation which shows that the CCR/CAR's are violated in the 3-point interacting model. 



II. A CLASS OF CHRONOLOGY VIOLATING MODELS 

In this section, we describe a class of nonlinear differential equations on "spacetimes" in 
which "space" consists of finitely many discrete points. By making identifications in these 
spacetimes, we introduce CTC's and obtain spacetime models generalising that studied by 
Politzer IJ. These identifications are implemented in the field theory by imposing certain 
boundary conditions which place constraints on the theory. 

Let S be a finiteQ collection of points S = {z a | a — 1, . . . , s} for some s > 2, and 
define spacetime to be the Cartesian product 5xR. We define fj to be the Hilbert space of 
complex- valued functions on S with inner product (/ | g) = J2 z es f( z )d( z )- This space has 
vectors v a as an orthonormal basis, where we define v a (zp) = 8 a p. With respect to this basis, 
we may write functions in fj as s-dimensional complex vectors, so that (/ | g) = pg = f a g a , 
where we sum over the repeated index. 

We will study model field theories derived from Lagrangians of form 

£ = I(V>ty - ^) - tfW*l> - ^W, (2.1) 

where ip(t) G $j, W is a self-adjoint positive operator on Sj and A G The corresponding 
field equation is 

^ = -iWijj - (2.2) 
which conserves the quantity ip^ip, and therefore reduces to the linear equation 

^ = -iWip-iX^ (0)V>(0)V>, (2.3) 
once the initial data ^(0) is specified. Thus the unique solution to Eq. (|2.2p with this data 



is 

rP(t) = e - iA ^ t(oWo) e-^V(0). (2.4) 

The configuration space variables ip a have conjugate momenta ^£,0 and Eq. ( |2.2| ) may 
be written in the Hamiltonian form 



1 Many of our results for free fields generalise easily to the case in which S is a manifold. 



2 Nai'vely, one might expect the momenta to be \iiplx- However, the Lagrangian ( |2.1| ) is a second 
class constrained system and the correct momenta may be obtained using Dirac brackets [20]. 
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with Hamiltonian 



h(if>, ifl) = £W a ^ p + ^W^pipa- (2.6) 

To introduce CTC's we partition S into two subsets Si and S 2 containing si and s 2 
elements respectively, with si + s 2 = s and s 2 < Sj, and make pointwise identifications of 
S 2 x {T+} with S 2 x {0 - } and S 2 x {T~} with 52 x {0+} for some T > 0. This idealises 
wormholes linking the lower surface of S2 at t — with the upper surface of 1S2 at t = T, 
and the upper surface of S 2 at t = with the lower surface of <S 2 at t = T. Note that CT 
and + (and correspondingly T~ and T + ) are regarded as distinct topological points for this 
purpose. 

The partition of S induces a partition of the basis vectors v a into the sets ei, . . . ,e si and 
fi, . . . , f S2 whose respective spans are denoted fji and fj 2 . Clearly, we have dimi^ < dimi^i. 
We will also write the projection of ip e onto fii {i = 1,2) as In the classical field 
theory, the identifications are implemented by the imposition of the boundary conditions 

MT-) =A^j 2 (0 + ) and ^ 2 (T+) = B^(0~), (2.7) 

where A, B are unitary maps of $) 2 to itself, corresponding to the evolution through the 
wormholes (cf. Goldwirth et al. f|0). Politzer fl[ takes A — B — II. The role of these 
boundary conditions is simply to ensure that the evolution round the wormholes is consistent. 
We require V ; i(^) to be everywhere continuous, although ipi(t) ma Y De discontinuous at 
t = 0, T. Thus (|2~2| ) is suspended at these points. 

Except in the special cases in which S consists either of 2 points or 3 points arranged in 
a ring, the interaction term in ( |2.1| ) is not a nearest neighbour interaction and is therefore 
rather unsatisfactory as a model field theory. We will therefore restrict our discussion of 
specific interacting models to these cases. A more realistic field theory will be discussed 



elsewhere 211. 



III. THE CLASSICAL INITIAL VALUE PROBLEM 

In this section, we examine the behaviour of the classical field equation ( |2.2| ) subject to 
the CTC boundary conditions (|2.7j ). For a generic class of W and T, we show that the free 
field initial value problem is well posed for data specified before the dischronal region. We 
then examine the nonlinear theory and show that (generically) solutions exist for all initial 
data specified before the dischronal region; moreover, this solution is unique in the case of 
"weak coupling" , but nonunique for "strong coupling" . 

To define the class of generic W and T, we decompose the operator e~ lWT (which imple- 
ments the free evolution between t = + and t = T~) in the block form 



<-''"'''=(££]■ (3-D 



3 In the notation of ||, our A is equal to (W^) 1 , and B is equal to W> . 
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with respect to the decomposition Sj = Sji © Sj 2 - The generic case is defined to be the case 
in which the normQ \\S\\ of S is strictly less than unity. (Note that \\S\\ < 1 because e~ lWT 
is unitary.) Equivalently, we require that Q should be an injection from Sj 2 into Sji so that 
Q has a left inverse K, (i.e., such that KQ = U|f, 2 ) which is uniquely specified if we require 
it to annihilate the orthogonal complement of ImQ. In this case we have \\K\\ = ||<5|| _1 - 
This requirement is the reason for our restriction that dim$j 2 < dim^h: otherwise, Q would 
necessarily have nontrivial kernel. The generic case corresponds to the situation expected in 
physically realistic field theories in which wavepackets spread out so that some proportion of 
any wave emerging from the wormhole at t = + manages to avoid reentering it at t — T~. 

A. Free Case 

We now show that, in the generic case with A = 0, the equation of motion (|2.2| ) with 
boundary conditions fl2.7| ) constitutes a well posed initial value problem for arbitrary data 
ip £ ft specified at t = (and therefore for any t < 0). In fact, we will only need the 
weaker condition that A — S be invertible on S) 2 - 

The evolution between t = + and t = T~ is given simply by the operator e~ lWT ; 
accordingly, given data at t — 0~, the problem reduces to the study of the evolution between 
t = 0~ and t = + . Because ip\{t) is required to be continuous at t = 0, it remains 
to determine ^2(0 + ) in terms of -^(O - ). The only constraint on ^(O" 1 ") is that the CTC 
boundary conditions be satisfied, i.e., that ip2{T~) = Ai/j 2 (0 + ). From Eq. ( |3.1|) we have 
ip 2 (T~) = i?-0i(O) + Sy>2(0 + ) so, provided A — S is invertible, ip2(0 + ) is uniquely specified as 

M0 + ) = (A-S)- 1 RM0). (3.2) 
For < t < T, the solution is thus 

^^- M 'i(A-s)%M0))' (3 ' 3) 

and in particular, we obtain 

^i(T)=M^i(0), (3.4) 

where the matrix M is 

M = P + Q(A — Sy 1 R. (3.5) 

The matrix M is unitary, as we now show. Let ij) 6 Sj± and use the unitarity of e~ lWT to 
compute 



4 If V\ and V 2 are complete normed linear spaces (Banach spaces) with norms || • ||i and || • H2 and A : 
V\ — > V2 is a linear mapping, then the norm \\A\\ of A is defined by \\A\\ = supj eVl ^ ^ ||A/1|2/||/||i- 



S 



UA-sy'R^- 



-iWT [ i> 
\{A — SY x m\) 

= \\(p + q(a - sy'Ryw 2 + 11(11 + s(a - sy^mw 2 

= \\Mij\\ 2 + \\A(A-Sy 1 R^\\ 2 . (3.6) 

By unitarity of A, we now have ||M?/>|| 2 = HV'II 2 and conclude that M is unitary. 

Thus we have shown that there is a unique classical solution for each choice of initial 
data at t = 0~ and that the full classical evolution from t = 0~ to t = T + is 



^(T H 



M 
B 



^>(o- 



(3.7) 



Moreover, the solution is clearly continuous in the initial data, so we conclude that this initial 
value problem is well posed for data specified in the past of the CTC region on surfaces of 
constant t. 

The situation is different for data specified between t = + and t = T~. Here, the 
initial value problem is well posed only for a subclass of data satisfying certain consistency 
requirements. For example, data specified at t = + must obey Eq. ( |3.2|) . This phenomenon 
has been noted before in various situations P,|5]J2^]; it arises because the CTC's introduce 
constraints on the dynamics and has important implications for the quantum theory. Note 
that one may nonetheless specify the data at any given point freely: it is always possible 
to choose the remaining initial data so as to satisfy the consistency requirements. Thus our 
system has a "benignity" property analogous to those discussed in P1J22|1. Related to this 
phenomenon is the fact that the classical evolution is nonunitary between t = 0~ and + 
and between t = T~ and T + . To see this, take any initial data with ipi(0) = 0: at t — 0~, 
the initial data has norm HV^O - )!!; for < t < T the solution vanishes identically; and 
finally, at t — T + , the solution again has norm ||V'2(0 _ )||. 

Finally, it is instructive to see how this classical evolution may be derived using the path 
integral methods of Goldwirth et al. . This is described in Appendix [A]. 



B. Interacting Case 

We now consider the full interacting classical field theory in the generic case. We will 
show: (i) there exists at least one solution for arbitrary initial data; (ii) there is a weak 
coupling regime in which there is a unique solution; and (iii) there is a strong coupling 
regime in which there exist many distinct solutions for each choice of initial data. 

In the absence of CTC boundary conditions, the solution is given by ( p.4[) . We write 
a = , 0i(O) and b = ^(O" 1 ") and implement the CTC boundary conditions by requiring b to 
satisfy 

Ab = e- iXT ^ a+m XRa + Sb), (3.8) 

for given a. 

To study the solutions to this equation, we first note that it implies ||6|| = \\Ra + Sb\\ 
and hence, by the unitarity of e~ lWT , that ||a|| = \\Pa + Qb\\. In the generic case (in which 
Q has left inverse K) any solution b must therefore take the form 
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b = K(U-P)a, 



(3.9) 



for some unitary U on Sj%. Substituting back into Eq. (|3Tq), and rearranging, we find that 
b solves ( |3.8|) if and only if 



KUa = Kf(U)a, (3.10) 
where f(U) = P + Q{A^V) _ S yi R and 

?7(I7) = ATa f {! + (£/- P)^K^K{U - P)} a. (3.11) 

Because t](U) is real-valued, /(£/") is a unitary operator on $j 2 . 

Clearly, any solution to the fixed point equation U = f(U) necessarily yields a solution 
to Eq. Q3.8p ; moreover, any such £7 must take the form U(z) = P + Q(zA — S)~ 1 R for some 
z on the unit circle. (Note that this expression is always well-defined because \\S\\ < 1.) 
Thus the problem of existence reduces to finding fixed points of the equation z = e lT! ^ u ^ 
on the unit circle. Now 

r](U(z)) = XT (\\a\\ 2 + \\(zA - S^Raf) , (3.12) 

which is a continuous single- valued function from the unit circle to the real line; thus e ir, ( u ( z " 
is a mapping of the unit circle to itself with vanishing Brouwer degree (see, e.g., Ch. 1 of [p3fl). 
Accordingly, for each choice of initial data a e Sj% there exists at least one fixed point of g 
and thus at least one solution to Eq. ([3.8|) , so we have proved the claim (i) above. 

To establish claim (ii), we write the right hand side of Eq. (|3.8|) as Ag(b) where g : S} 2 —* 
fj 2 and consider the fixed point problem b = g(b) on the ball B = {b G Sj 2 \ \\b\\ < r ||a||}, 
where r = ||<5|| _1 (1 + ||-P||). This ball contains all solutions to Eq. ( |3.8|) as a consequence 
of Eq. (|3.9|). For any b%, b 2 in £>, we have 

\\g{b x )-g{b 2 )\\ = ||(l-e aT (ll bl ll 2 ~ll fc2 " 2 ))(i?a + S6 1 ) 
+e iXT(\ M *-\ M *) s{bi _ b2) \ 

<( Co AT||a|| 2 +||5||)||6 1 -6 2 ||, (3.13) 

in which we have used the elementary estimates |1 — e ia \ < \a\ and | ||&i|| — H&2II I — \\bi — b 2 \\, 
and c = 2r (||i?|| + || S'H t^q) is a positive real constant depending only on P, Q, R and S. In 
the generic case, for AT||a|| 2 < c _1 (l — || || ) (i.e., weak coupling), g\ B is a strict contraction 
(which need not map B to itself) and standard contraction mapping arguments now imply 
that there can be at most one fixed point in B. Putting this together with (i) and using the 
fact that all fixed points of g must lie in B, we have proved (ii). 

Finally, to prove (iii) we note that if AT||i?a|| 2 ^> 1 (i.e., strong coupling) then the fixed 
point problem z = e %r,<yU<yZ ^ described above has many solutions on the unit circle; moreover, 
because Ra 7^ 0, these solutions must correspond to distinct values of b and hence of ipi(T). 

In Figs. [I] and |2|, we explicitly show how nonuniqueness arises in a model with two spatial 
points and P = —Q = R = S = l/y/2 with A = 1. For this model the classical solutions 
are in one-to-one correspondence with the solutions of the fixed point equation z = e ir, ( u ( z " 
on the unit circle, because K and U(z) are scalars. Writing z = e t6 , this becomes 
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((6) = 9 + 2kn, 



(3.14) 



for k G Z where = fi(l + \V2e ie — 1|~ 2 ) and we have written /x = AT|a| 2 forthe "coupling 
strength". This equation may be solved graphically by plotting each side separately for 
—7i < 9 < 7r looking for intersections. Fig. [I] shows the appropriate plots for fi = 0.5, from 
which it is clear that there is a unique solution, whilst Fig. |2| corresponds to the case /i = 3.0 
where there are 7 solutions. 



C. Preservation of the Symplectic Structure 

Except at t — 0, T, the classical dynamics is generated by the Hamiltonian h, and 
therefore preserves the symplectic structure on phase space in the initial and final chronal 
regions. Owing to the CTC boundary conditions, it is not clear that the evolution from 
initial to final chronal regions also preserves the symplectic structure. Here, we express the 
classical evolution in phase space language and prove that the evolution from initial to final 
chronal regions is implemented by a symplectic transformation. 

The classical phase space is T = C s with symplectic structure given by the 2-form 
Q = —idip^ A dip = —idip a A dip a . In the usual way, functions on Y are regarded as functions 
of independent variables ip and ip\ A symplectic transformation £ is a diffeomorphism of V 
which preserves Q, i.e., = Q, where is the pull-back of Q by £. This is equivalent 
(see e.g., §40 in p4| ) to the Dirac bracket relation 

{f^,goO D , x = {f,g} D , ${x) , (3.15) 

where the Dirac bracket of two functions on T is defined by 

st \ ( df dg dg df\ 

Corresponding to the decomposition ip = (ipx, ip 2 ), we have r = T 1 x T 2 , and associated 
natural projections tt^ : T —>■ T^. Then Q can be expressed as 

ft = TT^i + 7T*Q 2 , (3.17) 

where, for k = 1,2, Qk = —idijy k A dipt is the symplectic form on T^. Any unitary matrix 
U on i}2 defines a corresponding natural symplectic transformation of T%, which we denote 
Xu- in addition, for t G R, r t = exptldh is the evolution generated by the Hamiltonian h, 
where / is the canonical isomorphism between 1-forms and vector fields on T specified by 
ft(Jw, •) = oj(-). We have r 4 *fi = Q for all t. 

With these definitions, the diffeomorphism rj implementing evolution from t = 0~ to 
t = T + is rj = (k,xb), where B is the unitary matrix appearing in the CTC boundary 
conditions and k is a mapping from an open set U C Ti into Ti defined as follows. First, 
we define a differentiable map a : U -^Tasa solution to the equations 

TV i o cr = idi, 7r 2 o t t o a = xa ° 7T2 ° a i (3.18) 

which express the classical consistency requirement Eq. ( |3.8|) . Then k is defined by 
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K = Tli O Tx O 0~. 



(3.19) 



In general, there will be many possible choices for k reflecting the nonuniqueness of the 
classical evolution. For any such choice, the relation k*Qi = Q\ can be proved using the 
composition rule (/ o g)* = g*f* as 

K, ill = a TrpTT^li 

= (7*0 — a*T^7T2^l2 

= (7*«0i + tt*0 2 ) - a*n* 2 x* B ^2 

= (3.20) 

Thus, because x*b^ = ^2, we conclude that r\ = (k,Xb) preserves fl. 
In terms of Dirac brackets, writing ip(T + ) = r](ip,iip^), we have proved 

{MT + ),MT + )} D = 0, (3.21) 

and 

{VaT+),^(T + )} D = -«W. (3.22) 

Note that the evolution between t = 0~ and t = + (and similarly between t = T~ and 
t = T + ) is not symplectic in general. 



IV. THE QUANTUM INITIAL VALUE PROBLEM IN THE ABSENCE OF CTC'S 

In order to prepare for our discussion of chronology violating models, it is useful to show 
how a study of the QIVP reproduces the results of canonical quantisation for Eq. ( |2.2| ) in 
the absence of CTC's. We first discuss the case of Fermi statistics to avoid the operator 
domain technicalities of the bosonic case. 

The canonical approach starts by identifying the classical canonical coordinates ip a and 
iip^ and the classical Hamiltonian h^ipa^ipjj defined in Eq. (|2.6|). A Hilbert space $ is then 
constructed on which bounded operators ^1, . . . , \& s represent the CAR's for s degrees of 
freedom - that is, {^ a , ^> = and {^/q,,^/^} = 5 a p for all a, (3. The quantised (normal 
ordered) Hamiltonian H is defined as a (bounded) self-adjoint operator on £ by substituting 
^ a for ijj a in the RHS of Eq. (|2.6|) using its literal ordering. The quantum evolution generated 



by H evolves a general operator A from time to t by 

A(t) = e iHt Ae~ iH \ (4.1) 

and the evolved operator therefore satisfies the Heisenberg equation of motion 

A{t)=i[H ) A{t)}. (4.2) 

Thus, by virtue of the CAR's, ^ a {t) = e lHt ^ a e~ tHt solves the original equation of mo- 
tion ( |2.2| ) as an operator differential equation with initial data \I/ Q (0) = \l/ a . Moreover, the 
CAR's are necessarily preserved by this evolution. 



12 



It is possible to reproduce these results from a slightly different angle, namely by treating 
Eq. ( p. 2D as an operator differential equation and considering the quantum initial value 
problem (QIVP). Given initial data ty a representing the CAR's, we say that the QIVP is well 
posed if there exists a unique operator- valued solution ^f a (t) to Eq. (|2.2|) with \I/ Q (0) = *f> a 
and the evolution preserves the CAR's. To show that this is indeed the case, we note that 
for arbitrary initial data given as bounded operators on Eq. ( p. 2D has the unique solution 

^ a (t)=e- M ^{e- mt ) a ^p. (4.3) 

The proof of uniqueness closely parallels the analogous argument for the classical differential 
equation. One may check that this evolution preserves the CAR's either by explicit compu- 
tation or by noting that the above solution must agree (by uniqueness) with that obtained 
from the canonical approach. Thus the QIVP for Eq. (|2.2| ) is well posed in the fermionic 
case. 

Of course, it is not usually advantageous to consider the QIVP directly because it is 
rare that the equation of motion may be solved in closed form for general operator-valued 
initial data. However, for the chronology violating models considered in this paper, it will 
not always be possible to assume that the initial data is a representation of the canonical 
(anti)commutation relations and therefore the canonical method is no longer guaranteed to 
yield solutions to the equation of motion Eq. (|2.2|) . In these situations, we must therefore 
employ the more general setting of the QIVP. 

In the bosonic case, we encounter unbounded operators and therefore must proceed more 
carefully. We now describe the technicalities required in order to generalise the foregoing to 
this case.Q 

Definition Lei T> be dense in Hilbert space and let ^i(t), . . . , ^ s (t) be closed operator- 
valued functions on R such that V is a core for each ^ a (t) and is invariant under the ^ a (t) 
and^ a (t). Then the V a (t) are said to be a solution to Eq. j\2.Sj ) on V if each ^ a (t) is strongly 
differentiable with respect totonV with derivative — iW a ^ pit) — i\^ 1 (t)^ 1 (t)^ a (t) . 

Note that this definition extends that used above for the bounded case. 
Definition The closed operators . . . , \l/ s are said to represent the OCR's on $ if they 
have a common dense domain X invariant under both the ty a and the ty> a with 

= [*!,,*£]/ = (4.4) 

and 



5 For completeness, we give the following definitions (see pq] ). An operator A with dense domain 
D(A) in Hilbert space 5 is closed if its graph (i.e., the set of pairs (/, Af) as / runs through D(A)) 
is a closed subspace of J x J in the product topology; A is closable if the closure of its graph is 
itself the graph of an operator, called the closure of A. An algebraic subspace T> of J contained in 
D{A) is a core for a closed operator A if A is the closure of its restriction to T>. A densely defined 
operator A is essentially self-adjoint if its closure is self-adjoint. T> C D{A) is invariant under A 
if Af £ T> for all / G T>. Finally, an operator-valued function A(t) is strongly differentiable with 
respect to t on D with derivative B(t) if T> is contained in D(B(t)) and D(A(t)) for all r in some 
neighbourhood of t and Hfe -1 ^ + e) - A(t)) - B(t))f\\ -> as e -> for all / G V. 
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= <W (4-5) 

for all f G X and such that ^J^a (summing over the repeated index) is decomposable 
on X. That is, there exists a projection-valued measure Pq on R such that X contains 
V o = U^>o P [-^]^ and ^l^af = k^dP^f for all f EX. 



The reason for the technical requirement of decomposability is that it guarantees f26 
that all such representations of the CCR's are equivalent up to unitary equivalence and 
multiplicity (i.e., the conclusion of von Neumann's theorem holds). 

The canonical quantisation of Eq. (|2.2| ) proceeds as follows. Suppose that operators ^ a 
represent the CCR's on Hilbert space $ with dense invariant domain X, and let T>o C X 
be defined as above. The quantum Hamiltonian may be defined on T> by substituting the 
operators ty a into the RHS of Eq. ( |2.6| ) to yield an essentially self-adjoint operator whose 
closure is denoted by H. Moreover, Vq is easily seen to be invariant under e tHt for t 6 R. 
Thus, the evolved operators ^ a (t) defined by 

* a (t) = e im ^ a e- iHt (4.6) 

are strongly differentiate with respect to t on V with derivative ie lHt [H, ^/ a ]e~' lHt and the 
CCR's may then be used (on V ) to conclude that the ^ a (t) solve Eq. fl2.2|) on V Q in the 
sense defined above. 

As in the fermionic case, we may reproduce these results by studying the QIVP. The 
situation for general initial data is summarised by the following: 

Proposition Lei $ be a Hilbert space and V C $ be dense. Suppose that ty a , (a = 1, . . . , s) 

are closed (possibly unbounded) operators on $ such that 

(i) T> is a core for each *S? a and a domain of essential self-adjointness for 

(ii) T> is invariant under \£d and e~****i , *'r for all t e R. 

Then the operators ^ a (t) defined as the closure o/ e~ ,w 7*1 { e ~ lWt ) a 8 on ^ cons ^wte 
the unique solution to Eq. \2.%J on V with initial data 

An immediate corollary of this is that if the represent the CCR's on $ and the domain 
V is defined as above, then the QIVP for Eq. (|2.2| ) is well posed on V . 



V. THE QUANTUM INITIAL VALUE PROBLEM FOR CHRONOLOGY 

VIOLATING MODELS 



A. General Formalism 

We now analyse the quantum initial value problem for Eq. ( |2.2|) in the presence of 
CTC's, beginning with the case of the CAR's. Suppose that the operators ty a (a — 1, . . . , s) 
provide a representation of the CAR's on Hilbert space We specify these operators as the 
initial data for the QIVP at time t = 0~ . Writing and \&2,j to denote those operators 
associated with S\ and 52 respectively, we therefore seek operators \l/2j(0 + ) such that the 
evolution between t = + and T~ obeys the consistency requirement ^2,%(T~) = 2,j(0 + ) . 
Denoting vp^j = a;, \&2j(0 + ) = bj we therefore require the bi to satisfy 
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Remarkably, and in contrast to the situation for the classical theory, it turns out that this 
specifies the bi uniquely in the generic case as we now show. 

We first construct a solution to Eq. ( |5TT| ) and then prove its uniqueness. For z G C, 
let N(z)ij be the matrix- valued function of z defined by N(z) = (zA — S)^ 1 R, which is 
analytic in an open neighbourhood of the unit circle in the generic case. Then for any 
unitary operator V on Hilbert space we may use the (Dunford) functional calculus (see 
e.g., pp. 556-577 of |2"7fl ) to define N(V)ij as a matrix of bounded operators on &. Using 
this notation, Eq. ( |5.1| ) may be rewritten in the form 

b t = N(e iXT{a * ak+b * bk) ) ijaj . (5.2) 

Next, let $ r be the eigenspace of a\di with eigenvalue r and decompose # = ® r 3r- We 
emphasise that ajaj is not the total particle number on S at t = 0", but rather the particle 
number on S\. Thus, for example, £o is n °t 1-dimensional, but consists of all states at t = 0~ 
with no 5i-particles. We now define unitary operators U r on the 3V- by the recurrence relation 

U r+ i = exp i\Ta\ (s ik + N{U r %N{U r ) jk ^ a k , (5.3) 

with Uq = 1. Denoting U = ® r U r , it is easy to see that Eq. (|5.1| ) is solved by 

h = N(U) ijaj , (5.4) 



by comparing with Eq. ( |5.2| ) and using the fact that each aj maps $ r +i to 3v an d annihilates 

do- 

We now prove that (|5.4|) is the unique solution to Eq. ( |5.1|) . Suppose that b±, . . . , b s , 2 solve 



Eq. (|5.1|) , and write U = e lXT ( a k ak+b l bk \ Because N(U)ij is a matrix of bounded operators, 
Eq. ( |5.2| ) implies that the bi annihilate ^q. Accordingly, U leaves 5o invariant and U\$ = IL. 
Now suppose inductively that U leaves $ r invariant for some r > 0. Provided that r is not 
the largest eigenvalue of a\a k , Eq. ( |5.2| ) and its adjoint imply that each bi maps $ r +i to 
3V and each b\ maps 3> to $ r +i- Accordingly a\a k + b\b k and thus U leave 3y+i invariant. 
Hence by induction, we find that each ^ r is an invariant subspace for U, so we may write 
U = r U r with each U r unitary on $ r . It is then easy to see that the U r must satisfy the 
recurrence relation ( |5.3| ) with Uq = U\$ = 1. We have therefore completed the proof of 
uniqueness. 

Finally, we note that this solution is representation independent in the following sense. 
Suppose that ^ a form a Fock representation of the CAR's, and let 6, be the unique solution 
to Eq. (|5.1| ) on By the Jordan- Wigner theorem (see e.g., @), an arbitrary representation 
^>' a on J' takes the form ^ = U~ l {^! a ® 11) {7, where U : 5" ® 91 is unitary and 0T is 



an auxiliary Hilbert space. Then the unique solution to the analogue of Eq. (|5.1| ) on is 

6J = £/-!(&. ® !)[/. 

In the CCR case, certain domain questions must be addressed. We suppose that the 
ty a are a representation of the CCR's on $ with common invariant domain X and define 
Do C as in Section 0. An important property of this domain is that $ r C Vq for all r, 
where # r is again defined as the eigenspace of aja^ with eigenvalue r. Then it is easy to see 



15 



that the same construction as used in the CAR case yields a solution to Eq. (|5.1|) on V ; 
moreover, one may show that it is the unique solution such that V is a core for each 6 i5 and 
is independent of representation in the same sense as in the CAR case. 

Once the unique solution to Eq. (|5.1| ) has been obtained (for either CAR's or CCR's) we 
may substitute back to find 

and check to see whether or not this evolution preserves the CCR/CAR's and is therefore 
unitary. We will analyse various cases of this problem in the following subsections. 

B. Free Fields 

Here A = and Eq. (|5.2|) immediately yields the unique solution 

b t = (A-S)T l R ]k a k . (5.6) 
Substituting, we find that the evolution is given by 

*2,i(T + ) = Ai^(O-), (5.7) 

where M = P + Q(A — S)~ 1 R is unitary. Note that one obtains the same result for both 
Bose and Fermi statistics. This evolution is easily seen to preserve the CCR/CAR's; there 
is therefore a unitary X on $ such that 

# a (T+) = X^ a (0~)X. (5.8) 

An interesting feature of the above is that the operators hi are linearly dependent on the 
a;. Thus the components of ^(O 4- ) do not form a representation of the CCR/ CAR's for M+N 
degrees of freedom. In effect the system is reduced to only M degrees of freedom, reflecting 
the fact that the CTC's place N constraints on the system. Accordingly, the evolution 
between t — 0~ and t = + is nonunitary, although unitarity is restored at t — T + . In 
addition, we see that it is not legitimate to employ canonical methods to evolve the quantum 
field in the dischronal region (if one intends to solve the equation of motion Eq. fl2.2p ) because 
the data at t = + does not obey the CCR/CAR's. 

As a final check on our result in this case, and on the loss of degrees of freedom, let us 
quantise by the familiar method of obtaining classical mode solutions. Let ej(i) (respectively, 
fj(t)) be the classical solution to the free equation of motion with initial data ej(0~) = 
(fj(0~) = fj) where the basis vectors and fj were defined in Sect. [TT]. We write the 
quantum field ^(t) as 

*{t)=a i e i {t) + b j f j (t), (5.9) 

where the tij and bj form a representation of the CCR/CAR's on Hilbert space The 
components ^ a of the field are obtained by taking the inner product with v a . The time 
evolution of the a« and 6j is defined by re-expressing the field as 
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*(t) = a i (t)e i + b j (t)f j , 



(5.10) 



which leads quickly to the above unitary evolution from 0~ to T + using the results of 
Section JT|. In the dischronal region, however, fj(t) vanishes and so ^(t) = a^e^t) and the 
reduction to M degrees of freedom is explicit. 



C. Interacting Fields 

Here, we consider three simple examples. Model 1 is a system with two spatial points 
and yields a unitary theory for both Fermi and Bose statistics. Model 2 is a system with 
three spatial points. We study this theory for Fermi statistics and show that the resulting 
theory is nonunitary. For simplicity we work in the appropriate Fock representations and 
take A and B to be the identity. 

Model 1 Our set of spatial points is S = {z\, z 2 }, and Si = {^} for i — 1, 2. Thus W is a 
2x2 matrix and P,Q,R, S are scalars. 

Fermi statistics The Hilbert space $ for two fermionic degrees of freedom is isomorphic to 
C 2 <g> C 2 . The unique solution to Eq. ( |5.1| ) is 

b = (1 - S^Ra, (5.11) 

as is easily verified using the fact that e mata a = a. Substituting back, we obtain 

a(T) = (P + Q(l - S)~ 1 R) a, (5.12) 

which is identical to the unitary free evolution obtained in the previous subsection. This 
contrasts with the generically nonunitary evolution obtained by Politzer J^] for this model 
using the self-consistent path integral - see Section |IX. 

Bose Statistics Here, 5 = t 2 ®t 2 (where £ 2 is the Hilbert space of square summable sequences) 
and the unique solution to Eq. ( |5.1| ) takes the form 



b = f(a'a)a, (5.13) 

where / : N — > C is defined recursively by 

f(n + 1) = ( e ^r(»+i)(i+l/(«)l a ) _ sy'R, (5.14) 

with /(0) = (1 - S^R. 

Thus the evolution of a is given by 

a (T) = e -^at(l+|/(ata)P)a (p + Q /(a t fl) ) Q 

= e -iff(«ta) + Q( e *9(«ta) _ s^r) a, (5.15) 

where g is a real-valued function on N defined by g(0) = and g(n) = XTn(l + \f(n — 1)| 2 ) 
for n > 1. This may be rewritten as 

a(T)=X t aX, (5.16) 
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with X = e lh ( ata ) and h(n) denned by h(0) = and 

-ih(n+l) 



e -ih(n) e ~ig{n) 



X 



P + Q{ 



sy x R 



(5.17) 



The left hand side is always of unit modulus, so h(n) is real-valued and the operator X is 
unitary Thus the evolution from t = 0~ to t = T + is again unitary. We note that this 
theory agrees with the corresponding free theory on 5i (though the theories differ on $ r for 
r > 2). 

Model 2 In this example, our set of spatial points S = {z\, z 2 , z%}, is partitioned into 
Si = {zi, z 2 } and Si = {z%}. The matrix W is now a 3 x 3 self-adjoint, positive matrix, and 
the block decomposition of e~ lWT yields a 2 x 2 matrix P, a 2-dimensional column vector 
Q = (Qiy a 2-dimensional row vector R = (R 1 , R 2 ) and a scalar S. 

Fermi statistics The Fock space is 5 = C 2 <S> C 2 <S> C 2 . Given operators ai and a 2 at t = 0, 
we seek an operator b such that 



-iXT{a\ ai +a\a,2 +fe t 6) 



{Rid\ + R 2 a 2 + Sb) 



Using the results above, the unique solution to this equation is 



AXTa a 



S)~ 1 {R 1 a 1 + R 2 a 2 ) 



(5.18) 



(5.19) 



as may easily be checked by decomposing # = 5o © -Si © $2 with $ r the eigenspace of a\ak 
with eigenvalue r. 

Substituting, we find that 



01 (T) 
a 2 (T) 



a 1 

«2 



In Appendix 0, we show that 

(0 I (a 2 (T)ai(T) + oi(T)a 2 (T))o{aS | 0) = F (\Q 2 \ 2 - \Qi\ 2 ) detP, 



and 



where 



(0 I ai(T) 2 c44 I 0) = FQiQ 2 detP, 



F = e 
and a(S) is defined by 



AT 



i\T 



— s 



-i\a(S)T 
1 - S 



+ 



e -i\a(S)T _ 1 



1 



\s\ 



(5.20) 



(5.21) 



(5.22) 



(5.23) 



a(S) 



\Ri 


l 2 + 


P 2 


2 




l-5| 2 



11-51 



2 ■ 



(5.24) 



Thus, except in the free case or for very carefully tuned parameters the CAR's are neces- 
sarily violated and the evolution is therefore nonunitary. Note that the coefficients of F in 
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Eqs. ( |5.21 ) and (|5.22|) vanish simultaneously for all T if and only if W is block diagonal with 
respect to the decomposition i} = .fii © i} 2 (in which case \S\ = 1 and we are no longer in 
the generic case). 

In the next section, it will be useful to have an explicit matrix representation for this 
evolution. Choosing a basis for the 8-dimensional Fock space $ such that 



ai(0) 



10 \ 

1 







He 



a 2 (0) 



/ 



( 01 \ 



i o - 1 , 



is 



(5.25) 



where 1 2 is the 2x2 identity matrix, one may show that 



/ M l2 M a \ 

n Na 

U -N a 

J 



1U, 



(5.26) 



Here is the unitary matrix M = P + Q{1 — S) l R and N is also a 2 x 2 unitary matrix 
given by iV = (P + Q(e lXr — S)~ 1 R)U , where U is another 2x2 unitary defined by 



Un 7 & 7 



ij a j 1^2 



-i\T{a\a k +tfb) 



a i I $2 ' 



(5.27) 



(which makes sense because the exponential preserves $1). The precise form of iV will not 
concern us; however, we note that M 7^ N, because a\(T) and a 2 (T) fail to anticommute. 
Bose statistics The Fock space is £ 2 (g> £ 2 ® £ 2 and the unique solution to Eq. (|5.1| ) is 



b = f(dJd, c^c)c, 

where 

c = ||i?|| _1 (i?iai + i? 2 a 2 ), d = \\R\\~ 1 (R2ai — -Ri<z 2 ) 
and f(m,n) satisfies 

f(m, n + 1) = f e ^K(«+D(i+l/MP)l _ ^ _1 



\R\ 



(5.28) 



(5.29) 



(5.30) 



with f(m, 0) = ( f 
pendix |Cj that 



iXTm 



S) 1 \\R\\- Substituting back to determine Oj(T), we show in Ap- 



(0 I (a 2 (r)ai(T) - a!(r)a 2 (T))dt c t | ) = F (|Qx| 2 + |Q 2 | 2 ) det P, 



(5.31) 



with F given by Eq. (|5.23|) . This should be compared with Eq. (|5.21[) . Thus the evolution 
fails to be unitary on # 2 . 
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VI. DISCUSSION OF THE NONUNITARY EVOLUTION 



In the previous section, we showed that Model 2 was subject to a nonunitary evolution 
for both Bose and Fermi statistics. In this section, we discuss this evolution in more depth 
in the fermionic case. Recall that the Fock space 5 is 8-dimensional, and that the operators 
^a(0~) (a = 1,2,3) represent the CAR's for three degrees of freedom on Writing ^ 1)f 
for the operators associated with points in «Si, and ty 2 for the operator associated with 
the single element of S 2 , we write ^^(O - ) = tij for i = 1,2. The Heisenberg evolution 
*«(0") -> ^ a (T + ) is such that *i, f (T+) = a^T) and ^> 2 (T + ) = ^ 2 (T~). Our principal 
results in this section are, firstly, that the Heisenberg picture evolution cannot be expressed 
in either of the forms 

V a (T + ) = X-^ a {0-)X, (6.1) 

or 

V a (T + ) =X^ a (0-)X, (6.2) 

for some operator X on ^; secondly, that the Heisenberg picture evolution does not admit 
an equivalent Schrodinger picture description in terms of a superscattering operator. In 
addition, we will discuss the problem of extending the evolution from that of the ^ a (0 _ ) to 
arbitrary operators on 

Firstly, then, we show that the Heisenberg picture evolution cannot be expressed in either 
of the forms Eq. (|6.1| ) or ( |6.2|) . The form Eq. (|6.1| ) is clearly impossible because it would 
entail {ai(T) , a 2 (T)} = 0, and we may dispose of Eq. ( |6.2j ) as follows. The explicit form 
of the a,i(T) given above shows that any such operator X would necessarily preserve the 
subspaces $o,$i and #2 of moreover, because 

OiCOta = AfijOj-lft, (6.3) 

where M is unitary, we conclude that X\^ is unitary up to scale. Then it suffices to note 
that 

{ ai (T),a 2 (T)} I 11) = X\ ai XX^a 2 

+a 2 XX ] a 1 )X I 11) (6.4) 

which vanishes because X preserves #2 and X\^ is unitary up to scale. Accordingly, we 
cannot cast the evolution into either of the special forms Eq. (|6.1|) or (|6.2| ). 

Secondly, we show that the Heisenberg picture evolution cannot be described by a super- 
scattering operator. Recall that a superscattering operator on the state space of a (separable) 
Hilbert space $ is a linear mapping $ of the trace class operators T(3 r ) on $ such that if 
p G T($) is a positive operator^] of unit trace, then $p is also a positive element of T(#) 
with unit trace. Thus, $ is a linear mapping of density matrices to density matrices, which 
need not preserve purity. If a superscattering operator $ describes the Schrodinger picture 



6 For our purposes, a "positive operator" means one which is non-negative definite. 
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evolution of a system, then the Heisenberg picture evolution is given by the linear mapping 
$' of the bounded operators £(§") on defined by 



Trp($'Z) = Tr($p)Z, 



(6.5) 



for all p G T{$) and Z G £(#). In fact, $' is the dual mapping to $ under the natural 
identification of £(#) with the dual space of T(#). 

The dual mapping $' possesses three easily established properties: (i) $'1 = II; (ii) 
($'Zy = $'(Z j; ) for all Z\ and (iii) $' is positive in the sense that %' Z is a positive operator 
whenever Z is. If one writes the superscattering operator using the index notation (e.g., jl8] ) 



($P) A B 



then $' may be written as 



(%'Z 



I ry\D 



c 



v BC P Di 



BC 



(6.6) 



(6.7) 



Returning to our case of interest, we now show that there is no superscattering operator 
$ for which the Heisenberg evolution can be written as 



y a (T + ) = $'* a (0" 

We will need the fact that if u = (ui,M2) T and v 
column vectors and a G M, then the eigenvalues [i of the matrix 



(6.8) 

(^i,f2) T are 2-dimensional complex 



K 



( a U2 u\ \ 

«2 a v\ 

Mi a —V2 

\ uf —V2 a ) 



1: 



(6.9) 



satisfy 



(fi — a) 4 — (fx — ct) 2 (v)u + v'v) + \u'v 



t„|2 



0. 



(6.10) 



Note that the operator C = all + (3ai + f3a\ + 7a 2 + 'jal takes the form (|6.9| ) with u 
(P, 7) T and therefore has eigenvalues 



H = a ± VW+W- 



(6.11) 



Accordingly, C is positive if and only if a > \/\(3\ 2 + I7I 2 . Now suppose that there exists a 
superscattering operator $ such that Eq. Q6.81 ) holds. Then using the properties (i) and (ii) 
of $', we have 



S'C = at + !3a 1 {T) +/3a 1 (T) t 
+ 7 a 2 (T)+7a 2 (T) t , 



(6.12) 



which may be seen to take the form (|]|) with u = ((/3,7)M) T , v = ((J3, l)N) T . Hence its 
eigenvalues are 
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H = a ± vW+T)fTA, (6.13) 

where the two ± signs are independent and 

A = y/(\0\ 2 + N 2 ) 2 - |^| 2 (6.14) 

is real and positive by the Schwarz inequality and the unitarity of M and N . Because these 
matrices are unequal, we may choose (3 and 7 so that A > 0. Choosing a to lie in the range 

VW+\l\ 2 < « < vW + |7l 2 + A, (6.15) 

the operator C is positive, but $'C is not. Accordingly, $' violates property (hi) above and 
therefore cannot be the dual of a superscattering operator. 

Next, we consider the Heisenberg evolution itself in more detail. It is worth pointing out 
that we have not by any means obtained the full Heisenberg picture evolution; at present 
we know the evolution of only a 3-dimensional subspace (spanned by the \I/ Q (0 _ )) of the 
64-dimensional space £(#) of linear operators on the 8-dimensional Hilbert space Owing 
to our results above, various natural strategies for extending this evolution to the whole 
of £(5) are denied to us: the evolution cannot be extended as a *-homomorphism (i.e., 
mapping any polynomial in the \I/ a (0 - ) to the corresponding polynomial in the \E' Q (T + )) 
because the CAR's are violated; we cannot write Z — > X~ l ZX or Z — > X^ZX because of 
our observations above, nor can we write Z — > %'Z for some superscattering operator $. 

It therefore seems that there is no natural extension of our evolution to £(50- As a 
concrete illustration of this type of behaviour, let us consider an example with one fermionic 
degree of freedom. Define the operator a on C 2 by 



1 




(6.16) 



and suppose an evolution is given such that 11 — > II, a — > /xa and a* — ■> licl\ where < /x < 1. 
It turns out that there are at least two choices for the evolution of a^a consistent with 
a superscattering operator description. The first is that a^a — > a^a, corresponding to a 
superscattering operator $ with action 



a (3 \ fa ix (3 
(3 1 — a J I n/3 1 — a 



(6.17) 



on the state space of C 2 , whilst the second is a^a — > aa 1 " and corresponds to the superscat- 
tering operator £ with action 

W'T^Y (6-i8) 



(3 1 — a J \ ix(3 a 

To conclude this section, we note that the failure of positivity which showed the nonexis- 
tence of a superscattering operator can be traded for a loss of the trace preserving property: 
by allowing 11 — > kIL for k > y/2, any positive C of the form discussed above is mapped to 
a positive operator. One might therefore attempt to extend this in some way to a positive 
evolution on the whole of £(5) (which can be done if the evolution on 11, Oj(0), aj(0)t is 
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completely positive - see Theorem 1.2.3 in Arveson | |29|| ) thereby obtaining (by duality) a 
Schrodinger picture evolution possessing all the properties of a superscattering matrix ex- 
cept the preservation of trace. Rather than allowing individual probabilities to be negative 
with total probability equal to unity, we would now have positive probabilities with a total 
in excess of unity. It would be tempting to rescale this total to remove this problem, but 
that would amount to rescaling aj(T), for which there is no obvious justification. 



VII. THE CLASSICAL LIMIT 

With the normal ordering used above, we have shown that the quantum theory is uniquely 
determined in the generic case for all values of the coupling constant A. On the other hand, 
we have also seen that the classical theory is nonunique in the strong coupling regime. It is 
therefore interesting to determine the extent to which the classical theory may be regarded 
as a limit of the quantum theory. 

We consider Model 1 with Bose statistics. Reintroducing the units of action by replacing 
a and b by h l l 2 a and h l / 2 b respectively, the consistency requirement Eq. ( |5.1|) becomes 

b = e- ihXT ^ a+b]b) (Ra + Sb), (7.1) 

in which a and a) obey the CCR's [a, at] = 1. The unique solution to this is b = f(h\T] a' 1 a) a, 
where 

f(v- n +l) = ( e -(«+i)(i+l/(^.)| 2 ) _ sy l R, (7.2) 

with f{v] 0) = (1 - S)- 1 R for all v. 
One expects the classical limit to be 

b d = g(XT\a d \ 2 )a ch (7.3) 

where 

g(fi) = lim f(ji/n;n), (7.4) 

which implements the limit h —>■ while keeping XThn = \x fixed, where we define \i = 
\T\a c \\ 2 . Note that the classical quantities a c \ and b c \ have dimensions of (action) 1 / 2 . If the 
limit g(fj) exists uniformly in \x in some neighbourhood, then we have 

bd = (e iXT ^ 2+ ^ - S)' 1 R, (7.5) 

for initial data corresponding to this neighbourhood. This is, of course, equivalent to the 
classical consistency condition Eq. ( |3.8| ) for this model. 

Thus it is important to clarify the convergence properties of g n (n) = /(///n; n) as n — > oo. 
This turns out to be somewhat delicate and we have not as yet brought the general case 
fully under analytic control. However, we have investigated the problem numerically in the 
particular example P = —Q = R = S = l/y/2 and A — B — 1 for various values of /x. With 
these parameters, it is easy to see that (gn(fJ') 1 + 1) / "s/2 lies on the unit circle for all n, fj,; 
accordingly we define 
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(7.6) 



where — tt < arg z < n. In Table |, we give the values of 9 n for n = 500, 1000, 1500, 2000 and 
2500 with fj, = 0.5, 3.0, 7.0, and 13.0 and also give the value of 9 corresponding to the nearest 
solution to the classical equation of motion. These results suggest a slow convergence to a 
unique classical limit for each of these values of \i. At these values of /x, the classical system 
has 1, 7, 13 and 25 solutions respectively (see Figures [I] and |^ for the cases /i = 0.5 and 
fi = 3.0) of which the quantum theory picks out precisely one. However, it is by no means 
the case that 9 n converges for any value of /i. Indeed, it appears that convergence occurs 
only for certain bands of /i, the first two of which are < \i < 3.5 and 6.7 < \i < 8.5. We 
believe that this band structure continues indefinitely; it seems reasonable to assume that 
the bands become narrower as \x increases. 

One may interpret these results as indicating that for those ranges of the coupling 
strength \x = XT\a c \\ 2 in which the classical limit exists, the unique classical solution isolated 
by the quantum theory is the correct classical description of the dynamics. However, for 
other values of the coupling strength there is no classical limit and presumably no valid clas- 
sical description of the dynamics. Finally, we note that these results are strongly dependent 
on the choice of operator ordering as we will see in the next section. 



VIII. OPERATOR ORDERING 

So far, we have worked with a single choice of operator ordering, namely the literal 
ordering of Eq. ( |2.2| ) which corresponds to normal ordering of the quantised Hamiltonian. 
In this section, we briefly discuss the effect of allowing alternative orderings in which Eq. ( |2.2|) 
is ordered as 

jj, = -iWip - ia\(^ip)ip - i(l - a)\ip(ipty), (8.1) 

for a G [0, 1]. The foregoing treatment is the case a = 1. 

Consider Model 1 for Bose statistics. The analogue of Eq. (|Q| ) is 

b = e~ iaXT ^ a+b]b \Ra + S 6 ) e -*(l-<*)AT(ata+6t 6)) ^ 

(we have set A = 11 for simplicity). Making the ansatz b = /(a^a)a, we find that / satisfies 

/(0) = (R + Sf(0)) e-^Ci-^d+l/Co)! 2 ), (8.3) 

and 

f{n + 1) = e -^T(n+l)a(l+|/HP) ^ + Sf{jl + 

xe -iAT(n+2)(l-a)(l+|/(n+l)| 2 ) ; /g^ 

for n > 0. The case a = 1 was treated in Section and uniquely determines f(n + l) in 
terms of f(n) for each n. However, the case a = is rather different and is described by 

f(n) = e -^(™+ 1 )( 1 +l/(")l 2 ) (R + Sf(n)) , (8.5) 
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where we have written h explicitly. It is easy to recast this into the form of the classical 
consistency requirement Eq. ( p.8|) and it follows that f(n) is uniquely determined for small 
quantum numbers nh <C (AT) -1 but not for nh ^> (AT) -1 , i.e., classical nonuniqueness 
reemerges at high quantum numbers. There are therefore many functions f(n) solving 
Eq. ( |3.5| ), each one of which corresponds to a different "branch" of the quantum theory. 
Most of these branches do not possess a classical limit. However, in contrast to the situation 
for normal ordering, every classical solution will arise as the classical limit of some branch 
of the quantum theory. 

It would be interesting if the nonunitarity of Model 2 could be removed by a suitable 
ordering prescription. In Appendix |C], we investigate this for orderings of form ( |S.ip with 
the ansatz b = f(cftd, c [ c)c with c and d given by Eq. (|5.29| ). (Solutions do exist which take 



this form.) For all < a < 1, we find that the (anti) commutation relations are violated for 
generic values of the parameters. 

IX. SELF-CONSISTENT PATH INTEGRAL 

A. General Formalism 

In this section, we compare the results obtained from the QIVP with those obtained 
using the self-consistent path integral developed by Thorne and collaborators P,[7|,|50| and 
employed by Politzer |l| . To establish our notation, we briefly review the quantisation of our 
system by path integral methods in the absence of CTC's. Starting with the bosonic case, 
it is convenient to use the holomorphic representation (see, e.g., [|31|]) in which the Hilbert 



space # is the space of analytic functions f(c\, . . . , c s ) with inner product 

(/ | g) = j VciVce- ctc W)g(J) } (9.1) 

where we write to denote (cT, . . . , cj) and the measure is 

p c t Pc = TT^^. ( 9. 2) 

i 



The Hilbert space $ carries a (Fock) representation of the CCR's in which ct acts as 
multiplication by cj and Cj as d/dcj. Operators on $ are described by their kernels: 

(A/)(ct) = J Pc' t Pc'e- c ' tc 'A( C t;c')/(c' t ). (9.3) 

In particular, if K is a s x s matrix, then the mapping /(c^) — > f(c*K) has kernel expc^-fTc'. 

Starting with the (normal ordered) quantised bosonic Hamiltonian H on ^, one may 
obtain the kernel for U = e~ tHt in the form 

U t (J;c') = J n^7(t , ) t ^7(Oexp|i(7 t W7W+7 t (0)7(0)) + ^[ 7 ]}, (9.4) 
where the action functional 5 [7] is defined in terms of the classical Hamiltonian ( [2.6|) by 
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%] = I ( - 2 W)Hf) - - H(j(t')Mt'y) ) dt\ (9.5) 



and the paths ■jit') are subject to the boundary conditions j'(t) = and 7(0) = d. In the 
free case, for example, one may evaluate the path integral explicitly to give 

U t (J;d) =expc t e-* m c'. (9.6) 



One may develop the path integral treatment for Fermi statistics in a parallel fashion |3T| 



by replacing the integration variables by Grassmann numbers and VcDc^ by Berezin mea- 
sure. Again, the resulting kernel has the action of e~ lHt on J, where H is now the fermionic 
normal ordered quantised Hamiltonian. 

A natural generalisation of this to enable the treatment of chronology violating systems is 
the self- consistent path integral [J|ffl,[7],M. Instead of integrating over all field configurations 
with 7(0) = d and 7^(T) = c' to form the kernel Uj<(c', d), the self-consistent path integral 
prescription requires that one should restrict the class of field configurations to those obey- 
ing the self-consistency requirements imposed by any CTC's present (here, the boundary 



conditions (|2.7|) ). To implement this, we first decompose $ = #1 <S>#2? where #1 is the space 
of analytic functions in variables CL\ , . . . , Cl Sl , and $2 is the space of analytic functions in 
b\, . . . , b S2 . The (self-consistent) evolution kernel from t = 0~ to t = T + can then be written 
in the form 

X{a\ &t ; a ', b') = Ue h]Bb 'UT{a)- a'). (9.7) 

Here, AT is a normalisationconstant and the factor e b]Bb ' implements the boundary condition 
ip2{T + ) = Bip2(0 ) while Ut is given by the same path integral as Ut but taken over all field 
configurations with 7(0) = (a',b"), J'(T) = (a,Ab"y for any b". As noted by Politzer [|TJ, 
Ut(cl'', a!) may be obtained from Ut{o^ , b'; a', b') by setting b = Ab' and integrating over all 
possibilities in the Hilbert space measure of # 2) that is, 

U T {o)-a') = J VtfVbe- btb U T (a^tfA^;a',b), (9.8) 

which may be rewritten in the form 

%{a^\al) = J VtfVb J VdVce- bU - cic 

xe btA ' c U T (a\J;a',b). (9.9) 

Thus, by expanding e 6tAtc as 

(9,0) 

we obtain the matrix element (m | Ut | m') in the form 

(m I U T I m') = 

(m; n | Ut | m'; n), (9-H) 
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where the vector | m) e $ 1 is the function n«( m iO 1 ' 2 W ni i and the vector | n) G #2 is the 

function UiM'^iAtyi* ■ In addition, (A*b)i denotes YiiAjh- We refer to E 1- ( BD > 
which is a generalisation of the expression given by Politzer |]] as the partial trace definition 
of the self-consistent path integral. 

The fermionic case follows a similar pattern, when one replaces the integration variables 
by Grassmann numbers and uses Berezin measure; the main difference lies in the partial 
trace definition. Starting from the analogue of fl9.9|), we expand e fotAtc as 



e^=n(i+4(^c); 

8=1 

-E(- 1 ) n Il( Atc )" i (&Ir> (9-i2) 

n i 

where n = Yli n ii an d therefore obtain 

(m I % I m') = ^(-l) n (m; n\U T \ m'; n), (9.13) 

under the assumption that the Grassmann number b\ commutes with the kernel of 
Ut{cl \ c<; a', b), which holds if H conserves particle number (as it does in our case of in- 
terest). The factor of (— l) n was omitted by Politzer 0; it arises because terms of the form 
b\{A^c)i coming from e btAtc must be rearranged in order to move the fet's into the ket and 
the (A'c)i's into the bra of the matrix element (m; n | Ut \ m'; n). In Appendix [FJ, we will 
see how, for free fields, these factors ensure that the evolution computed from ( |9.13|) agrees 



with that obtained directly from the path integral, and also with that obtained from the 
QIVP. 



B. Free Fields 

Whilst one can use the partial trace definition to compute the quantum evolution X for 
free fields (see Appendix |B]), it is easier to evaluate the path integral directly, using the fact 
that the kernel of the free evolution is given by 

[/ T ( c t ; c ') = exp Je- iWT c'. (9.14) 
Writing e~ tWT in the block form (|3.1f) as above, we obtain 

U T {a ] ] a) = J Vb ] Vb exp {-b\t - A ] S)b + a) Pa' 

+a)Qb + tfA ] Ra'} , (9.15) 

which may be evaluated to give 

Ur{a)]a') = (det(l — A^S))^ 1 exp a 1 Ma', (9.16) 

where M = (P + Q(A — S)~ 1 R). In the generic case, the convergence of the path integral 
is guaranteed because H^S 1 !! < 1 and so 1 — A* S has positive hermitian part. 
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Noting that V(a)\ a') = expa^Ma' is the unitary kernel, because M is unitary, we con- 
clude that the unitary kernel obtained from the self-consistent path integral is 

X(a\ 6 + ; a', b') = exp {a^Ma' + tfBb'} , (9.17) 

whose corresponding operator X acts on annihilation operators and 6j according to 

X ] a i X = M ij a jl X%X = B^bj. (9.18) 

Moreover, the normalisation constant is given by M = det(l — A'S). 

In the fermionic case, the path integral may be evaluated explicitly to obtain a unitary 
evolution with the action ( |9.18|) on annihilation operators and normalisation constant M = 
det(l - 

Thus in both cases, we have obtained agreement with the QIVP evolution. Moreover, 
we have given a general proof of the unitarity of free field evolution using the self-consistent 
path integral; previously this had only been established in a particular case [l|]. 



C. An Interacting Model 

We study Model 1 of Section |V| for both Bose and Fermi statistics, employing the partial 
trace definition, and choosing the normalisation constant so that (0; | X | 0; 0) = 1, which 
is reasonable because the Hamiltonian H is particle-number preserving. In the fermionic 
case, we obtain 

(0 | Ut I 0} = (00 | e~ iHT | 00) - (01 | e~ lHT | 01) 
= 1-5 

(1 | Ut I 1) = (10 | e~ iHT | 10) - (11 | e' lHT | 11) 

= P - (PS - RQ)e~ iXT , (9.19) 

from which it follows that the evolution from t = 0~ to t = T + is given by 

(m; n \ X \ m; n) = 5 nn >5 mm >f(m), (9.20) 



where 



1 m = 

f{ m ) = \ P-( PS -RQ)e-^ T _ i (9-21) 



l-S 



m 



Thus X is nonunitary in general, which is essentially the result obtained by Politzer |I[ in 
special cases, modulo some changes of sign owing to the factors of (—1)™ discussed above. 
Except when \T/(2tt) G Z this differs from the unitary evolution obtained from the QIVP. 
In the bosonic case, we have 



x> e ~i\T(m+n)(m+n-l)/2 



(m\U T \ rri) = Y . , ni/9 (00 | (Ra + Sb) n (Pa + Qb) m | m'n), (9.22) 

z — ' [mini ) L / Z 

n=0 v 1 

and therefore conclude that (run \ X \ m'n') = 5 mm >8 nn r f(m) with /(0) = 1 and 
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(9.23) 



One may show that X fails to be unitary in general. Again, it clearly differs from the unitary 
evolution obtained from the QIVP. 



X. CONCLUSION 

In this paper, we have analysed in detail the classical and quantum behaviour of a class 
of nonlinear chronology violating systems. Classically, we found that unique solutions exist 
for all choices of initial data in the linear and weak-coupling regimes, whilst the solutions 
become nonunique in the strong-coupling regime. This confirms the expectation that the 
behaviour of nonlinear fields interpolates between that of classical linear fields and hard- 
sphere mechanics. Quantum mechanically, we have shown that one can make sense of the 
quantum initial value problem for chronology violating systems; moreover, (at least with a 
natural choice of operator ordering) the quantum dynamics is unique for all values of the 
coupling constant. We have also exhibited examples in which this evolution does not preserve 
the (anti) commutation relations; it seems highly likely that this is the general situation. 
Moreover, the nonunitary evolution cannot be described by a superscattering operator - the 
loss of unitarity is more radical than previously thought, e.g., by Hawking |fL8 |. 



We have also compared our quantum evolution with that computed using the self- 
consistent path integral, and found that they do not agree. This is not surprising, be- 
cause the equivalence of these approaches for non-chronology violating systems relies on the 
existence of a foliation by Cauchy surfaces and there is no a priori reason to expect the 
equivalence to persist in the presence of CTC's. In this regard it is interesting that the 
QIVP and self-consistent path integral are nonetheless equivalent for linear fields. To some 
extent, it is a matter of taste which approach one prefers. We prefer the QIVP approach 
for two main reasons. Firstly, we have found circumstances (e.g., Model 1 in Section |V C| ) 
in which one obtains a unitary theory from the QIVP but not from the path integral. Sec- 
ondly, the effect of the CTC's in our models is to introduce constraints which lead to a 
nontrivial geometric structure in the classical phase space. We suspect that the quantisa- 
tion of this system requires more than just a restriction of the class of allowed histories, 
and that the path integral measure should also be modified (a similar comment has also 
been made in Jl6|]). A hint of this appears in the treatment of linear fields, in which the 



propagator obtained from the self-consistent path integral must be rescaled by a factor of 
det(U — A^S)^ 1 . It is plausible that in the linear case, the required modification to the path 
integral measure reduces to rescaling by this constant factor, but that for the nonlinear case 
the modification is nontrivial. At present it is not clear to us exactly how the path integral 
should be modified; on the other hand it is clear that the QIVP does correctly implement 
the CTC constraints and remains close to the spirit of the classical treatment. 

The relationship between the unique quantum theory and the nonunique classical the- 
ory is intriguing. We have seen that there exist ranges of the coupling strength in which 
the quantum theory has a classical limit which selects precisely one of the many classical 
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solutions, and other ranges in which no classical limit exists. It would be interesting to 
understand the underlying reasons for this behaviour in more detail. 

Finally, it is curious that the classical symplectic structure can be preserved for systems 
which do not preserve the quantum commutation relations. It is tempting to wonder whether 
there is a way of quantising these models so that unitarity is preserved. Our uniqueness 
result for the QIVP rules this out within a Hilbert space context (at least with normal 
operator ordering) but it is possible that the situation might be different for the QIVP on 
an indefinite (Krein) inner product space in which irreducible non-Fock representations of 
the CCR's exist for even a single degree of freedom |J2 . The motivation for studying Krein 



spaces would be that the loss of physical degrees of freedom in the nonchronal region might 
be equivalent to the addition of unphysical states with negative norm-squared. 
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APPENDIX A: PATH INTEGRAL APPROACH TO THE FREE CLASSICAL 

EVOLUTION 



In this Appendix, we show how the classical evolution derived in Section [111 A| may be 



reproduced using a method due to Goldwirth et al. || and based on path integrals. (Gold- 
wirth et al. regarded the classical wave equation as the first quantisation of an underlying 
particle mechanics.) The central idea is to sum the propagators for all possible trajectories 
through the CTC region. We will use this method to determine the propagator between 
t = 0~ and t = T + , essentially repeating the calculation of || in our (slightly simpler) 
notation. 

The block matrix decomposition Eq. ( |3.1| ) suggests that we break the problem into four 
parts, evaluating the propagators from <Sj at t = CT to Sj at t = T + separately for each 
i, j = 1, 2. Note that a particle on 5 2 at t = 0~ must enter the wormhole there and reemerge 
on S2 at t = T + . Thus the S2 — ► S2 propagator equals B, whilst that for S2 —>■ Si vanishes. 
In addition, the propagator S\ — > S2 also vanishes by the time reverse of this argument. It 
remains to compute the propagator for S\ — > S%. In this case, there are countably many 
possible trajectories. The particle can either go directly to «Si with propagator P, or it can 
enter the CTC region to arrive at S2 at t = T~ (propagator R), pass through the wormhole 
to 1S2 at t = + (propagator A' 1 ), execute n circuits of the CTC's (propagator (A^ 1 S) n ) and 
finally travel from S2 at t = + to «Si at t = T + (propagator Q). The combined propagator 
for this trajectory is Q(A~ 1 S) n A~ 1 R; summing over all possible winding numbers and the 
direct trajectory, we obtain the total propagator 
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m = p + qI ^(a- 1 ^)™ a^r 



v n=0 



P + Q(A - S^R, (Al) 



which agrees with the result obtained in Section [III A 



APPENDIX B: PARTIAL TRACE FORMALISM FOR FREE FIELDS 

In this Appendix, we derive the evolution operator for free field models in the presence 
of CTC's using the partial trace formulation of the self-consistent path integral. 

We consider a general free theory whose Fock space is built using creation operators 
a[, . . . ,a\ and b\, . . . ,b\ , acting on vacuum | 0; 0). The aj and hi obey the CCR/CAR's. 
The basis elements are writtenF] 

| m; n) = \\[m^l\a\)^ f[(n t \)-^ (b\r I 0; 0), (Bl) 



i=l i=l 



and we write m = J^m,, n = J2 n i e ^ c - We will also need the alternative basis | m; n) 
defined by 



Sl S2 



m; n) = nK ! r 1/2 («!r n^ ! )- 1/2 ((^)lr I 0; 0>. (B2) 



i=l i=l 

Suppose the evolution U on Fock space is unitary and such that 



cii(T) = U ] aiU = Pijdj + Qijbj 

bi(T) = U j kU = Rijttj + Sijbj, (B3) 



where the matrix 



is unitary. We note that this entails that U preserves the total particle number a|aj + 

We now specialise to the bosonic case. From Section [IX A| , the evolution operator X has 
matrix elements given by 

(m; n | X \ m'; n') = AT b 5 n n' 5^(m; n" | U \ m ; n"), (B5) 



7 We note in passing that the basis used in Eq. (2) of Ref. Q for fermionic systems is not properly 
anticommuting. 



31 



where A/b is a normalisation constant, chosen to ensure that (0 | X | 0) = 1 (as it should be 
for any free theory). This allows us to evaluate N\> explicitly, because the matrix element 
(0;n | U | 0;n) is 

(0;fi| f /|0;n) = (0;0|n^^|0;n) = (0;0| f /n i 4^^ ( B6 ) 

and is therefore equal to the coefficient of \\ i X T i n the expansion of Yli(J2j(^ ^)ij x j) ni ■ 
We have used the fact that U preserves the vacuum. The generating function G(x\, . . . , x S2 ) 
for these coefficients can be found in §66 of |33[|, and is given by 



{-l) S2 (xix 2 . ..x S2 ) :I 
det(AtS' — diag (x± , x^ , . . . , xjj-)) 



(■lr.. • • • , x S2 ) = ^ , 7 —. (B7) 



The sum over all n of these matrix elements is obtained simply by evaluating the generating 
function with all X{ equal to unity. Thus we obtain 

A/; = det(H - A*S). (B8) 

Next, we claim that 

X^aiX = M ijaj , (B9) 

where M = P+Q(A — S)~ 1 R is unitary. Together with the trivial evolution X~ 1 b i X = Bijbj, 
this shows that X is unitary. Moreover, this is the free evolution derived in various ways in 
the body of the paper. 

To establish (|B9D, we first note that 



^^(m; n | Ubi | m'; n) = ^^(m; n | (A*b)iU | m'; n) 

n n 

= ^(m;n | UA^R^ + S jk b k ) \ m'; n), (BIO) 

n 

where the first step follows by relabelling the sum over n*. Collecting terms in the bi and 
rearranging, we have 

^(m; n\Ubi\ m'; n) = ^(m; n | U(A - S)^Rjka h | m'; n), (BH) 

n n 

and hence 

y~](m; n | aiU | m'; n) = ^(m; n \ l 'il\,<i, + Qijbj) | m'; n) 

n n 

= ^(m; n \ UM^aj | m'; n), (B12) 

n 

where M = P + Q(A — S)~ 1 R. Thus we have a{X = XM^aj as required. 
In the fermionic case, we define the operator X by 
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(m; n | X | m'; n') = Af<W ^(-l) n "(m; h" \ U \ m'; n"), (B13) 

n" 

where Af is chosen to ensure that (0 | X | 0) = 1. The factor of (—1)™" is necessary in order 
to obtain agreement with the canonical theory. To see this, note that the first step in ( |B10|) 
is not valid in the fermionic case, due to the anticommutation relations satisfied by the aj 
and hi and the definition (|B1|). Instead, the corresponding result is 



J2(-l) nim+m,) (m; h\Ubi\ m'; n) = ^(-l) n(mW) (m; n | (A%U | m'; n), (B14) 

n n 

in which the factors of (— l) m and (— l) m ' arise from anticommuting bi past the string of 
creation operators for | m) and | m') respectively. We may replace (_i) n ( m + m ) by (— l) n 
because U preserves the total particle number and therefore the summands can be nonzero 
only when m' = m + 1. 

Exactly analogous arguments to those for the bosonic case then show that eq. (|B9|) holds, 
and that X is unitary. Thus we have obtained agreement with the canonical theory. 

The constant Af is easily evaluated once it has been expressed in the form 

Af 1 = (0; | [/\ S2 (U - A^S)b S2 . . . 6 X ] b\ . . . b\ 2 | 0; 0), (B15) 
for then one may use the exterior algebra definition of the determinant^] to conclude that 

Mi= [det(l- A^S)]'\ (B16) 
To establish Eq. ( |B15| ), we write its RHS as A/" -1 and expand the exterior power to obtain 

AT- 1 = ^(0;0 | (-l)"^^...^ 1 ^!...^ | 0;0), (B17) 

n 

where c[ n ^ is defined to be equal to bi if rij = or (A^Sb)i if rij = 1. Next, move the leftmost 
cf 1 ^ with iii = rightwards using the anticommutation relations until it sits next to b\, at 
which point the bib\ combination may be removed by a further application of the CAR's. 
Repeating the process until all c^' s have been removed, one eventually finds 

A"- 1 = J](-l) n (0; | {A*Sb)Y | 0; n), (B18) 

n 

which is easily shown to be equal to ^2 n (—^) n (0;n \ U \ 0; n) = J\ff , thus verifying our 
claim. 



8 Recall that if Q is an endomorphism of a vector space with anticommuting basis bi, . . . , 6jv 5 then 
(A N Q)b N •••&! = (Qb) N ■ ■ ■ (Qb)i = (det Q)b N . . . b x . 



33 



APPENDIX C: VIOLATION OF CCR/CAR'S IN THE 3-POINT MODEL 

We present here the details of the calculation leading to Eqs. ( |5.21| ), ( |5.22| ) and 



and the statements made at the end of Sect. VIII. We consider the 1-parameter family of 



operator orderings labelled by a G [0, 1] discussed in Sect. |VIII| for which 



(t) = e -iAT a vf(0)t*(0) [ e -iW $ | J e -iAT(l-a)vK(0)t*(0) ; (q-Q 

and consider Bose and Fermi statistics simultaneously, seeking solutions of the form b = 
f($d, $c)c where 

C = \Tn\\ all d " = j-j-z— jj (C2) 

||ii|| ||ix|| 

obey the same commutation relations as and a 2 - 

Applying b to elements of form (c^) m c^ | 0), we obtain the consistency requirement 

f(m,0) = (\\R\\ +5/(m,0))e- iAT[m+(1 - a)(1+l/(m ' 0)|2)] , (C3) 

and applying 6 to elements of form (S) m (c^) n+2 | 0) for m,n > 0, we obtain the recursion 
relation 

/(m,n+l) = + S'/Crrz,/! + 1)) e -^[ m +o,(™+i)(i+|/(^)| 2 )+(i-o ; )(«+2)(i+| / (^+i)P )] _ (C4) 
We compute the quantity 

J = (0 | (a 2 (l>i(T) =F ^(T^T))^ | 0) (C5) 
for Bose (— ) and Fermi (+) statistics and also 

K=(0\ ai(T)W | 0) (C6) 

for Fermi statistics. Note that d)c^ | 0) is an element of $2; in the fermionic case it is 
a\a 2 I 0). 

First note that 

(0 I a t (T)a,(T)Sc^ | 0) = e^(0 | (P ik a k + Q^e-^^^^Pjm + Qjb)d^ | 0), (C7) 
where u = XT [(1 -a) (2+ 0)| 2 ) + 1]. We have 

P u Ri , P n R 2 -Pi 2 Ri, 
\\R\\ \\R\\ 

- ^- * detP - d , (C8) 



11*11 S||fl|| 
where we have used the identities 

PijRj = -QiS, (C9) 
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and 



P n R 2 - PnRl = -9^L, (cio) 

which follow from the unitarity of e~ lWT . The second of these is derived by noting that 
SRj = —QtPij. Similarly, we have 

Q 2 S Q 1 detP 

P l 2i<H = —T l — ^CH — = a, Cll) 

||i?|| S\\R\\ K } 

in which we have used P21-R2 ~ -P22-R1 = Qi det P/S. 
We can now compute 

,-,,detP f . ,„ „ / S 



(0 I a 2 (TK(T)rft c t I 0) = e"^^ {i^ 2 (/(1,0) 



\R\ 



|Q 2 | 2 (/(0,0)-^e-™°)l 2 }, (C12) 

where the ± is + for Bose and — for Fermi statistics. To obtain J, we interchange the 
suffices 1 and 2 and add for Bose and subtract for Fermi. The sign reversal occurs because 
d)c* I 0) flips sign under the interchange. The final result is then 

f f(l 0) - fro nW iAT l/(°i°)l 2 e -iAT|/(o,o)| 2 
J = e- det P(±\ Ql \* + |Q 2 | 2 ) I f J±°> + £ - — - } . (C13) 

The matrix element K may be computed for Fermi statistics as 

( f(\ Q) _ f(Q n V- iAT l/(0,0)| 2 -iAT|/(0,0)| 2 _ 1 

K = e-^Q.Q, det P \ f[ ' ' f T ' + -— } . (C14) 

\ S\\R\\ \\R\\ 2 ( 

With the particular operator ordering used in Section |V C| (a = 1), we have /(0,0) = 
117211(1 - S)- 1 and /(1,0) = \\R\\{e iXT - S)~ x . Substituting these values into Eqs. ( PT3|) 
and (PH ) and using the fact that \\Rf = 1 - \S\ 2 , we obtain Eqs. ([BTMl) , flpp and EH) 
respectively. 

Finally, one should also check that the expression enclosed within braces in Eqs. ( |C13|) 
and ( |U14 ) does not vanish. For AT <C 1, one may prove this by perturbing about the free 
solution to obtain /(0, 0) and f(l, 0) to second order in AT if S ^ R. If S is real, one needs 
to go to third order. 
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TABLES 



TABLE I. Table of 9 n defined by Eq. (7.6) for different values of the coupling strength \x. The 
final row shows the value of 9 corresponding to the unique classical limit picked out by the quantum 
theory. 



n 


H = 0.5 


H = 3.0 


6„ 


\i = 7.0 


H = 13.0 


500 


0.9020 


-2.7496 




2.2226 


2.7558 


1000 


0.9018 


-2.7492 




2.2124 


2.7527 


1500 


0.9017 


-2.7490 




2.2115 


2.7516 


2000 


0.9016 


-2.7489 




2.2112 


2.7510 


2500 


0.9016 


-2.7489 




2.2109 


2.7507 


Classical solution 


0.9015 


-2.7487 




2.2101 


2.7494 
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